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Abstract 

II ' A method for calculating the pion structure function directly in 

Mh. terms of light-cone wave functions is suggested. Taking twist-2 and 

r~ I ! twist-4 pion light-cone wave functions into account, it is shown that 

the QCD sum rule prediction is in agreement with quark distribution 

obtained from analysis of the Drell-Yan process. Twist-4 quark-gluon 

light-cone wave functions give a large positive contribution to the pion 

C^ ' structure function for xb < 0.2. A new constraint on the twist-2 pion 

light-cone wave function is obtained. We argue that the leading twist 
pion wave function f^iu) ~ 1 for u = 0.3 with an accuracy of about 
20-30%. 
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1 Introduction 

The calculation of structure functions for deep-inelastic lepton-hadron scat- 
tering, which is equivalent to finding the quark and gluon distributions in 
hadrons, is one of the most important problems of quantum chromodynamics 
(QCD). Perturbative treatments of deep-inelastic scattering give important 
information about structure functions as a function on Q^ [jl[ (Q^ is the mo- 
mentum squared of the virtual photon) and, under some assumptions, can 
produce reggion behavior at small xb = Q^ I^Vi) HI- At the present time 
there is significant activity in the application of the perturbative expansion 
for structure functions 0]. These perturbative investigations are very impor- 
tant. However, since quark and gluon distributions can not be determined 
from perturbative calculations, these distributions have been taken from ex- 
periment. 

Theoretical calculations in QCD have achieved some progress in determin- 
ing the second moments of the quark distribution functions for the nucleon 
[|] and pion |^ from QCD sum rules. The first moment of the chiral-odd 
structure function h\ was evaluated in |^. Recently, higher moments of 
quark distribution functions have also been considered in the QCD sum rule 
approach |]^. 

In ref . P] , loffe proposed a four-point correlator for theoretical calculations 
of quark distribution functions in the QCD sum rule approach suggested by 
Shifman, Vainshtein and Zakharov |^ . This method was applied for calcula- 



tions of nucleon structure functions, such as Fi{xb) [^j 9i{^b) and (72 (xg 



[|TT|], hi{xB) HTU- In the present paper we use a light-cone QCD sum rule 



[|T^. Recently, it was suggested to apply this version of QCD sum rules to 
the problem of calculating the deep-inelastic structure function [l^. This 
light-cone QCD sum rule is based on the fact (which was noted by loffe in 
[§]) that if Xb is not close to the boundary xb = and a;^ = 1, then the 
imaginary part of the deep-inelastic scattering amplitude is determined by 
small distances in the t-channel. In the case of the pion structure function, 
the nearest singularity in the t-channel for the correlator of two vector cur- 
rents and one axial current with a pion in the initial state is at t = — Yzf~P^ 
for highly virtual photons, where p is a momentum of the axial current. This 
means that in the case of intermediate xb we can apply operator product 
expansion (OPE) for the calculation of the correlator in the Euclidean re- 
gion. Then, we can use a dispersion relation to construct QCD sum rules. 



These light-cone QCD sum rules are formulated in terms of so-called light- 
cone wave functions of hadrons introduced in perturbative QCD to describe 
hadron form factors at large Q^ [|T3], [IB], |r^ . 



This new version of QCD sum rules can be a useful tool for investigations 
of quark and gluon distributions and light-cone wave functions in hadron 
physics. Furthermore, in principle, it is possible to incorporate perturbative 
corrections to improve the calculations. 

In the present paper we do not take into consideration perturbative cor- 
rections. So, our results can be considered as an initial condition for the 
evolution equation |]1|] . For this reason, at the end of this paper we compare 
our results with the quark distribution function for low Q^ |]l8l obtained from 
the the analysis of the Drell-Yan process. 

We consider the limit of massless quarks iriq = 0. In this limit, the pion 
is also massless, m^ = 0. 



2 Correlator 



As in [141, we consider the correlator 



T,,x{p.q^) = ~i j d'xd^ze^^^^''^^ < Q\T [j^x) , j%z) , jtmTi' {k) > (1) 

for our calculation of the pion structure function. Here k is the pion momen- 
tum, 

Jl = M7/.75C^, Jp,A = dlp,xd, (2) 

and the following kinematics is used: 

e = 0- q^ = {p + q-kf- t={p-k)^ = 0; s = {p + qf- 

Q^ = -q^; {2k,p + q) = s + Q^; {2pk) = p\ (3) 

The discontinuity in s at fixed p^ and Q^ of the correlator (|I|) is calculated 
from 



ImTnox = — 



Tf,px{p ,q ,s + ie)- T^px{p ,q ,s-ie) , (4) 



where p"^ and q^ are space-like vectors, p^ < 0, q^ < 0, such that |p^|, |g^| S> 
^QCD- In the scaling limit, we assume that \p^\ <^ \q^\ and keep only the 
first nonvanishing terms in an expansion in powers oip^/q^. 



ImT^pX is calculated in the physical region of the s-channel, and the pion 
contribution in this amplitude has the following form, 



/mT^pA = pj-^lm [zjd'ze^^^ < 7r{p)\T{ji{z),fM}Hk) >} • (5) 



On the other hand, the general form for the imaginary part in (|^) is 



Iml^tj d'ze''^^ <7r{p)\T{ji{z),jtm\7f{k) > 

+A,is, Q')iq'g,, - g«g« - gf gf + g«gf ) 
+Bi{s,Q)\p —Q^'j \P —1j 

^n( n^^( P^^'^ (i)W P'J^"^ m\ 

+B.(..0^)(p-^,'^')Jp-^,'^')^. (6) 

where q^^^ = q and g'-^-* = p + q — k are the momenta of the virtual photons; 

(g(l))2 = (g(2))2 = q2_ 

The structure function of deep-inelastic scattering is defined when q^^^ = 
q^'^\ and it is clear in this case that 



Q 

TT 

where 



Axlq\xB) = — (B,is, g2) + B^is, Q^) + B;{s, Q') + B^{s, Q')) , (7) 

TV \ J 0^-+00 



xlq'^ixB) ( pq \ ( pq \ . .gx 

Here q'^{xB) is d-quark distribution function of a pion. 



In this paper the tensor structure p^PpPx in correlator (0) is considered. 
We define the imaginary part of the correlation function for these tensor 
structures as fT,-^^t{p'^,XB)- Then, the dispersion relation for the function 
t has the following form, 

where the last term corresponds to the higher-states contribution. 

To suppress the contribution of exited states, as usually done in QCD 
sum rules, we will consider instead of t{p'^,x) its Borel transform in p^, 

f 2\n+l / ^ \ " 

t{M',XB)= , lim /~^ -^ tip^xs) 

'xIq^^b) + / p(s, XB)e-'l^''ds^ . (10) 

The left-hand side will be calculated in terms of the light-cone pion wave 
function by using the operator product expansion (OPE). Every new term 
of the OPE will be suppressed by a factor {Aqcd/MY^'^, where t is the 
twist of light-cone wave function. For the right-hand side, we will use the 
standard continuum model for higher states, whose contribution is suppressed 
exponentially. As usual in the QCD sum rule approach, we can find the 
desired physical quantity (here it will be the quark distribution function) by 
using a fitting procedure in the region for the parameter M ^ where the higher 
states (which are estimated in the continuum model) and higher-twist pion 
wave function contributions are small. 



3 Twist- 2 light- cone wave function. 

In this section we consider the contribution of the leading twist-2 pion light- 
cone wave function in the QCD sum rule ([T0|) . 

It is clear that in the formal limit when \p'^\, Q^ -^ oo, we can use a free 
(i-quark propagator. All interactions with soft nonperturbative gluon fields 
are suppressed by a factor A.qcd/p^ ^^ ^qcd/Q"^- Perturbative contributions 
should be taken into account. However, as pointed out in the introduction, 
we do not take these corrections into consideration, assuming that our results 



can be compared directly to the quark distribution function at low Q^, where 
these contributions are small. 

The result of a very simple calculation with free d-quark propagators in 
(HD gives 



T,,,{p,q,k)=tld''xd^ze'P^^"^'- 



(x — z)aZi3 



Att^{x - zYz^ 

< 0\uixht^l5lalplt3lxdiO)\Tc{k) > . (11) 

According to eqs.(^,|^), the interesting contribution corresponds to the sym- 
metric part of the amplitude T^pA- This symmetric part of the product of 
7-matrices in (|TTp has the following form, 

7M757a7p7/37A = ^p/3(^pa7A + 
9\a'y^^) + g(3xigf,alp + gpalp) 

+ {terms with gpx,gp^„ gx^,). (12) 

Clearly the terms with gp\,gp,_i,g\^ do not contribute to the Bi that are 
defined in eq.(|3). 



From eq.(|T2|) it is easy to see that the amplitude ( |TT|) is determined by 



only one matrix element, 

nl 
/O 



< 0|u(x)7;.75C^(0)|7r(A;) >= ikpf^ j dve-'^'^^^i^^iv) + x^g^{v) + 0{x')) 

2 \ rl 



+/- (^M-lr^J 1^ due-'^^'^^ig^iv) + 0{x'')),{l?,) 



where (pT^{y), gi{v) and g2{v) are twist-2 and twist-4 light-cone pion wave 
functions. In the formal limit \p'^\ — >■ cxd the leading contribution is deter- 
mined by the twist-2 pion wave function, ip^i'^j)- All other terms in the 
expansion of the amplitude (0) are suppressed by (Aq(-,^/p^)*^^ (t is the 
twist of a light-cone wave function) . 

Using the definition (p!3[), we obtain the following expression for the con- 
tribution of the leading twist-2 wave function. 



U t dv I d'^xdh '^^(^) e^(p-fc^)^+^(p+'?-^-^)^ 
■'Jo J in^x^z^ 

[zpiXf^kx + xxk^) + zx{x^kp + Xpkp)] . (14) 
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Here we have made the shift x -^ x+z. After integration over the coordinates 
X and z we obtain 



-f'lj^ 



ipTr{v)dv 



{p — kvY[p + g — kvY 
[{p + q- kv)p{{p - kv)f,kx + (p - kv)xkf,) 
+ (p + g - kv)x{{p - kv)^kp + (p - kv)pk^))] . (15) 

Using the following relations, which hold for the kinematics @, 

{p-kvf =p^{l-v)] {p + q-kvf = s- {s + Q^)v, (16) 

we obtain an expression for the coefficient of the tensor p^^PpPx in the basis 
(p, q^^\ g*^^)) (see ||), namely, 

Jo p^[s — (s + (y^'j-y) 

From eq.(|l^, it is easy to find the imaginary part of this amplitude, tijP', xb), 



as 



tjP^^XB) _ {l-v)(p.„{v) 

_ V9^(l -Xb) 
Q^p^ ■ 

In eq. ([l8D we have used the relationships 



(18) 



Comparing the result (|T8|) with the dispersion relation (|]), (and using 
isotopic symmetry, (pTt{v) = ip^il — v)), one may be tempted to claim that 

q{xB) = <^t,{xb)- (20) 

Relation ( [20| ) corresponds to a pure parton picture when the pion consists of 
two free quarks, since 

/ q{x)dx = 1; / xq{x)dx = 0.5, (21) 

Jo Jo 



which follow from normalization (/q (pT^{v)dv = 1) and the symmetry y:>Tr{v) = 
(fTr{l — v) of the twist-2 pion wave function. Note, however, the identification 
( POD can not actually be valid, since the light-cone wave function is not a 
positive definite function. From the point of view of the QCD sum rules that 



are being constructed in this paper, (|T9D means that the leading-twist pion 



wave function determines the quark distribution function only in the region 
where (Pt,{xb) has a positive value. 

Note that in the parton model, the quark distribution function is equal to 
(p'^{xb)- Thus, we can expect that the minimal corrections to the relation (pO|) 
(due to the higher-twist light-cone wave functions) will occur in the region 
where ipn{u) ~ fHu) ~ q{u) ~ 1. From the quark distribution function 
obtained in |TB[, one can find that q{xB) — 1 for 0.2 < xb < 0.3 with an 
accuracy of about 20%. So, we can expect that the QCD sum rules will be 
the most reliable in that region. To check this assumption, we will have to 
evaluate the contribution of twist-4 light-cone wave functions. 

4 Twist-4 Wave Functions 

The easiest part of the calculation is the contribution of the two-particle 
twist-4 wave functions gi and g2 (see (|13|)). The calculations with these wave 
functions are technically the same as those presented in the previous section. 
This correction was calculated in WM , and the corresponding contribution to 



t{p'^,XB) has the following form, 

w 2 N 4 ( gi{xB) + G2ixB) 1 dgiixBY 

U{p ,xb) = — + -g2{xB) 



p^ \ Xb 2 dx 



B 



= -,U^b)- (22) 

There are additional corrections corresponding to three-particle twist-4 
wave functions. The general form of these three-particle wave functions can 
be written as follows: 

< 0\u{x)gsGai3l'rl5{z)d{0)\n{k) > 

= Ukagpy - kpg^,)jVaie-'^-^^-''''^-'f{a,) 



(kx) 






k^{z^kp - zpK) Jva,e-''^'''-''''''Ua,), (23) 



i < 0\u{x)gsGai3{z)'y^d{0)\'7i(k) > 
-(^k^i^akp - xpk^)lva,e-'^^'''-''''''-^f, 



OLi 



k,{z^kf, - zpk^) J Va,e-''^''^-'''^'''Ua,). (24) 

Here and below we have used the notation JT>ai = /q daida2da3S{l — ai — 
^2— «3)- Comparing with the standard form for the hght-cone wave functions 
(when z = vx) that were considered in ||T9|, we obtain the following relations, 

f{ai) = ip±{ai), 
fxi^i) + fz{oii) = Vwic^i) + V5±(«i), 

f{ai) = ip±{ai), 

fxioci) + fz{ai) = <^||(ai) + <^±(ai)- (25) 

For practical calculations, it is very convenient to use the following for- 
mula for the quark propagator in the presence of the gluon field |jl3[ , 

<0\T{q{x)l,ql{z)}\0>=6" 



2'K^X^ 

+t^''^'''^^f^duG%{ux + {l-u)z) 

((1 -U){x- 5)7^7j, + U7^7,.(X - Z))af3 



ix — zY 



(26) 



X — z 



-g&^^' / du{x - z),Al{ux + (1 - u)z)—^ + 0(/). (27) 

JO ^ ZTT^lx — zr 



The contribution of the term (pq) to the diagram where a d-quark interacts 
with gluon field between points x and z can be written in the following form, 

3271^ Jo J {x-zYz^ 

< 0\u{x)'Jf,'j5Go,f3iuX + (1 - U)z){{l -U){X- Z)'^alp 

+M7a7/3(£ - z))7p^7AC?(0)|7r(A;) > . (28) 

9 



The result of our calculations with (pS]) is 
w 2n _ ^^ T If a2du'Dai 



— 1 /"i f^-^B da^ 

p^ Jo Jo xb + [i- — uja^ 



-{2(p±{ai) + (fll{ai)]a2=XB,ai=l-a2-a3 



p4 



• (29) 



In this equation we have used the fact that our result depends on the functions 
/x5 /z, fx and fz in the combination (|25|) only. 

The contribution of the diagram in which a quark interacts with soft 
gluon fields in the propagator < q{z)q{^) > is suppressed by a factor 1/Q^. 

The contribution of the last term (pTf) of the quark propagator is 



1 
dv{x - z)^A^{xv + (1 - v)z)d{Q)\'n{k) > . (30) 



In the Fock-Schwinger gauge {x^A^{x) = 0), the gauge field A^ can be ex- 
pressed in terms of G^^j as 

Ai,(xv + (1 — v)z) = / udu(vx + (1 ~ v)z)sGsuiu(xv + (I — v)z)). (31) 
Jo 

The result of our calculations for (ISO) is 






In eq.(P^) we did not include the part that could not be represented in terms 
of the light-cone wave functions ip\\, ip±, ip\\, (p^. This contribution comes 
from the following integral, 

^ ^^ x^z^ [ik,x + z)^ ik,vx + z))- ^ ^ 
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Numerical calculations show that the contribution of twist-4 quark-gluon 
wave functions (without the contribution of (p3D) cancels the contribution of 
the two-particle wave functions in the QCD sum rule for the first moment 
of quark distribution function when we use a self-consitent set of twist-4 
light-cone wave functions (they will be defined in the next Section) with the 
parameter e = 0. It can be proved that there are no higher-twist corrections 
to the sum rule for the first moment of the quark distribution function. This 
indicates that the term (^) should be small, and we can not exclude the 
possibility that its contribution is equal to zero for this particular set of 
light-cone wave functions. In any case, we do not have a model for the wave 
functions fx and fz that we can use to make realistic numerical estimates for 
the contribution (^). 

The remainder of the contribution (|30D has the following form, 

2 /-i-^'s das 

Jgiy^B) — ~ I '^\\\(^i)\ax=a2-a3;a2=XB 

p JO tt3 

~7a ^^ ¥'||(ai)f^ai|a2=i-ai-a3 (34) 

p JO 0:3 Jl—XB—013 

5 Light-Cone QCD Sum Rule. 

Now we can begin our analysis using the QCD sum rule for the quark distri- 
bution function (p^O]). This QCD sum rule has the following form, 

^n{xB) - Th(M^b) + fgixs) + fgl{xB)) = q'^ixs) + c{x B)e~'^''^''^' . (35) 



Here we use the results (p!8| , |22| , |29| , |34|) obtained in the previous sections. The 
last term on the right-hand side of (^) imitates the higher-states contri- 
bution, where m is the mass of a resonance. Below, we will assume that 
m ~ rriA^ = 1.25GeV. To find the quark distribution function q{xB) we have 
to analyze the QCD sum rules for all values of xb- 

In Fig.l, the dependence of functions {fi{xB) + fg{xB) + fgii^B)), fii^XB) 
and {fg{xB) + fgi{xB)) on xb is shown. The QCD sum rule is reliable in the 
region 0.2 < xb < 0.6, where the contribution of the twist-4 wave functions 
is small. Recall from our earlier discussion that we expected the minimal 
contribution of higher-twist wave functions to occur when 0.2 < xb < 0.3. 
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>< 



f4('<B)+fg(''B)+f.,(XB) 

f4(''B) 

fE(XB)+tg,(XB) 



Figure 1: Twist-4 contributions to the d-quark distribution function as given 
in Eq. (35), evaluated with the pion wave functions given in Eq. (36). 



For our numerical estimates, we have used the following set for light-cone 
wave functions: 



ipM = 6m(1 -u), 

5 1 6 

g^(u) = -5W^ + -e6^[uu{2 + 13mm) + 10m^(2 - 3m + -u^)ln{u) 

6 

+10m^(2 - 3m + -M^)/n(M)], 
5 

92{u) = — 5^mm(m-m), 

G,{u) = ^^^mV, 
o 

(/9x(«^) = 306\ai - a2)al[^ + 2e(l - 2a3)], 

(p\\{a.i) = 1206'^e{ai - Q;2)Q;ia;2Q;3, 

^^{ai) = 305^al{l - ag)!^ + 2e(l - 2^3)], 
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V9||(ai) = -1206^aia2a3[- + e(l - Sas)]. (36) 

We have used the notation u = {1 — u). The wave functions of twist-4 are 
very numerous. Here we use the twist-4 wave functions with leading and 
next-to-leading conformal spin [0 (see also [pO|] ). 

One of the parameters in (|36D is defined by the matrix element 



< 'K\gsdGa^i'jau\0 >= i6'^ f-nq^,- (37) 



The QCD sum rule estimate of ref.[|Tl yields 5^ = Q^GeV"^ at /i = IGeV 



(/i is the renormalization scale). The last parameter is associated with the 
deviation of twist-4 wave functions from their asymptotic form. Here we used 
the set with e = 0. It was noted in the previous Section that we can expect 
that the contribution of ( P^ ) is small or even equal to zero. To use the set 
with e = 0.5 (see Ref. [jl9[) we have to take into account the term (^3)). 

Note that there is large uncertainty in the choice for the form of these 
light-cone wave functions. Even the twist-2 pion wave function ipT^{x) is not 
known well enough, and there are different models for it. The most peculiar is 
the light-cone wave function with a two-humped profile, which was suggested 
by Chernyak and Zhitnitsky ||15 . 



The asymptotic form for the wave functions can be taken as a starting 
point for a detailed investigation of these QCD sum rules. The results thus 
obtained from the sum rule ( pSj ) are shown in Fig.2. A comparison of these 
results with the "experimental data", which were obtained from Drell-Yan 



process and extrapolated to the low normalization point [^], shows good 
agreement. Agreement with the "experimental data" can be improved by 
using a different form for the light-cone wave functions. 

It is interesting to estimate the second moment of the quark distribution 
functions. Assuming that the region near the end points xb = 0, 1 (where 
our considerations are not valid) gives a small contribution to the second 



moment, we obtain the following sum rule from eq.(35). 



where 



1 5"^ 

- - -— + 0{M'^) = M^ + {higher resonances), (38) 



M^= j\Bq\xB)dxB. (39) 

Jo 
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QCD sum rule predictions 

quarlc distribution from the Drell-Yan process 

QCD sum rule prediction with f4 contribution only 

q(XB)=6xB(l-XB) 




Figure 2: Theoretical results for the d-quark distribution function of the pion 
compared to experimental results. The solid curve is the complete result 
appearing in Eq. (35), evaluated with the wave functions given in Eq. (36). 
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Note that the contribution of the quark-gluon wave function is small. It is 
about 10% of the contribution of twist-4 two particle wave functions. This 



means that our QCD sum rule is almost the same as the one obtained in W^ 
The numerical evaluation of this sum rule gives Mg = 0.27±0.05, which is 
in good agreement with "experimental data": M| ~ 0.3 (see ||T8||). 

It was noted above that from the point of view of the partonic model 
we can expect minimal higher-twist corrections to occur in the region where 
q{xB) — 1- The "experimental data" presented in Fig.2 show that for xb — 
0.3, q{xB) = 1 • The twist-4 correction is small for xb = 0.3. Thus, we 
can conclude that our consideration is self-consistent and leads to a new 
constraint on the pion twist-2 light-cone wave function, namely, 

ip^{u) ~ 1 {u = 0.3). (40) 



This constraint follows from the "experimental data" [|T8[ and from the sum 
rule (|35|). The estimate of higher-twist effects confirms our earlier argument 
in favor of such a constraint. 

At the present time there are three different models for the twist-2 pion 
wave functions. The first is asymptotic wave function, 

^l'-{u) = 6n(l - u) <^""(0.3) = 1.26. (41) 



The second is Chernyak-Zhitnitsky wave function ||15|| , 

(f^^iu) = 30u{l -u){l- 2uf; <^^^(0.3) = 1.01. (42) 

And, the third was suggested by Braun and Filyanov |]T9| , 

if^^{x) = 6a;(l - x){l + a2(/i)3/2(5(2x - 1)^ - 1) 
+a4(/i)15/8(21(2x - 1)^ - 14(2x - 1)^ + 1)) 

(^f^(0.3) = 0.72 (43) 

The coefficients 02 — 0.41, a^ = 0.23 correspond to the normalization point 
/i = l.SGeV (see [g^). 

According to our findings, the most preferable model for the pion wave 



function is asymptotic one (^H). 
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6 Conclusion 

In this paper we have presented a new version of the hght-cone QCD sum rule 
that can be used for investigating the structure functions of deep-inelastic 
scattering. We demonstrated that this approach gives predictions for the 
quark distribution function that are in reasonable agreement with present 
"experimental data". 

We noted that the present QCD sum rule provides the possibility to 
determine the numerical value of the light-cone wave functions in the region 
0.2 < X < 0.6, where it is expected that the QCD sum rule should work. 
However, we have to emphasize that the perturbative corrections may give a 
significant contribution. These perturbative corrections can be evaluated in 
formalism of string operators p^ . 



Our analysis using the sum rule ( p5[ ) can give important information 
about pion wave functions that are currently widely used in calculations of 



different hadronic processes p3| . A new constraint for the twist-2 light cone 
wave functions was formulated from this sum rule. This constraint indicates 
that the twist-2 pion light-cone wave function is not very different from its 
asymptotic form. 

The agreement with "experimental data" can be significantly improved 
by changing the parameter set for the the light-cone wave functions. It is 
important to evaluate perturbative corrections. 
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